x_max = 2*pi;
N = round(2*pi/0.01);
mu1 = -1;
mu2 = -10;
max_iter1 = 30;
max_iter2 = 5;
FTCS_max1 = zeros(1,max_iter1+1);
u = initial_data(N,x_max);
FTCS_max1(1) = norm(u,inf);
for i = 1:max_iter1
    u = FTCS_onestep(u,N,mu1);
    FTCS_max1(i+1) = norm(u,inf);
end

FTCS_max2 = zeros(1,max_iter2+1);
u = initial_data(N,x_max);
FTCS_max2(1) = norm(u,inf);
for i = 1:max_iter2
    u = FTCS_onestep(u,N,mu2);
    FTCS_max2(i+1) = norm(u,inf);
end
plot_fig1(FTCS_max1,FTCS_max2);